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ABSTRACT 

We present analytical modelling of conical relativistic jets, in order to evaluate the role 
of the jet opening angle on certain key parameters that are inferred from VLBI radio 
observations of blazar nuclear jets. The key parameters evaluated are the orientation 
angle (i.e., the viewing angle) of the jet and the apparent speed and Doppler factor of 
the radio knots on parsec scales. Quantitative comparisons are made of the influence of 
the jet opening angle on the above parameters of the radio knots, as would be estimated 
for two widely discussed variants of relativistic nuclear jets, namely, those having 
uniform bulk speed and those in which the bulk Lorentz factor of the flow decreases 
with distance from the jet axis (a 'spine-sheath' flow). Our analysis shows that for both 
types of jet velocity distributions the expectation value of the jet orientation angle at 
first falls dramatically with increases in the (central) jet Lorentz factor, but it levels 
off at a fraction of the opening angle for extremely relativistic jets. We also find that 
the effective values of the apparent speeds and Doppler factors of the knots always 
decline substantially with increasing jet opening angle, but that this effect is strongest 
for ultra-relativistic jets with uniform bulk speed. We suggest that the paucity of 
highly superluminal parsec-scale radio components in TeV blazars can be understood 
if their jets are highly relativistic and, being intrinsically weaker, somewhat less well 
collimated, in comparison to the jets in other blazars. 

Key words: blazars: general - galaxies: active - galaxies: jets - galaxies: nuclei - 
galaxies: radio continuum 



1 INTRODUCTION 

Although Very Long Baseline Interferometry (VLBI) moni- 
toring of the radio knots in blazar jets has revealed several 
sources containing knots with apparent speeds v apP in ex- 
cess of 25c (e.g., Kellermann et al. 2004; Piner et al. 2006), 
the typical speeds for blazars known to emit the highest 
energy 7-ray photons (TeV blazars) are found to be much 
more modest, with v apP < 5c, and their radio knots are often 
subluminal (e.g., Edwards & Piner 2002; Piner & Edwards 
2004). The glaring contrast with the very large bulk Lorentz 
factors in the parsec-scale blazar jets (L > 30c), as inferred 
from TeV flux variations (e.g., Krawczynski, Coppi & Aha- 
ronian 2002), has been highlighted by several authors (e.g., 
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Giroletti et al. 2004; Gopal-Krishna, Dhurde & Wiita 2004; 
Piner & Edwards 2004, 2005; Levinson 2006). Similar ultra- 
relativistic bulk Lorentz factors have also been inferred from 
the intraday radio variability of some blazars (e.g., Rickett, 
Kedziora-Chudczer & Jauncey 2002; Macquart & de Bruyn 
2006) and, more directly, from VLBI measurements of the 
brightness temperatures of several blazar nuclei (e.g., Ho- 
riuchi et al. 2004; cf., Kovalev et al. 2005). The large mis- 
match between the estimates of P derived from these differ- 
ent types of observations has led some authors to postulate 
a dramatic jet deceleration between sub-parsec and parsec 
scales (e.g., Georganopoulos & Kazanas 2003; cf., Piner & 
Edwards 2005). An alternative approach has been to invoke 
a 'spine-sheath' configuration for the jets such that the fast 
spine close to the jet axis is surrounded by an appropri- 
ately slow moving sheath (e.g., Baan 1980; Komissarov 1990; 
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Laing 1993; Meier 2003; Ghisellini, Tavecchio & Chiaberge 
2005). 

These considerations recently led us to undertake an 
analytical study which showed that the modest apparent 
speeds of the knots of blazars, which are mostly unresolved 
by VLBI, can be reconciled with the extremely relativis- 
ts bulk motion inferred for TeV, and some other, blazars 
as noted above, if one considers a modest full opening an- 
gle (u> ~ 5°-10°) for the parsec-scale jets (Gopal-Krishna, 
Dhurde & Wiita 2004, hereafter Paper I). We also showed 
that the actual viewing angles, 9, of such conical jets from 
the line-of-sight can be substantially larger than those com- 
monly inferred (e.g. Jorstad et al. 2005) by combining the 
flux variability and the VLBI proper motion data (Gopal- 
Krishna, Wiita & Dhurde 2006, Paper II). Direct support for 
the assumption of conical parsec-scale jets comes from the 
VLBI imaging of the nuclear jets in the nearest two radio 
galaxies, M87 (Biretta, Junor & Livio 2002) and Centaurus 
A (Horiuchi et al. 2006). 

Our focus in this paper is on making a quantitative com- 
parison of the predictions for two widely discussed jet forms 
(i.e., those with uniform T and those with velocities decreas- 
ing away from the jet axis, the spine-sheath types), of the 
extents to which some key parameters of the VLBI radio 
knots in blazar jets would be influenced by the jet opening 
angle. The three main jet parameters examined are the ap- 
parent speed (flapp) of the knot, its Doppler factor (8) and 
the viewing angle (9) to the jet axis from our line-of-sight. 
Although discussed specifically for blazar jets, the present 
results are also relevant for 7-ray burst sources (GRBs) 
which too are believed to arise from extremely relativistic 
jets, with P > 100, and with even larger solid angles (e.g., 
Meszaros 2002; Liang et al. 2006, and references therein). 



2 COMPUTATION OF THE APPARENT 
PARAMETERS FOR CONICAL JETS 

Following the analytical prescription of Papers I & II, we 
approximate a radio knot with a circular disc-shaped region 
of uniform intrinsic synchrotron emissivity, and thus uniform 
surface brightness, presumably arising from a shock in the 
jet (e.g., Marscher & Gear 1985; Hughes, Aller & Aller 1985). 
See Fig. 1. 

We define the angle between the axis of the jet and 
the line-of-sight from the stationary core of the blazar to 
the observer to be 9 and the velocity parameter of the knot 
along the jet to be (3 = v/c. Because we will consider jets 
with finite opening angles, and thus with a range of viewing 
angles to different portions of them, we use the more general 
form of the usual relations, 



8 = 



r(l-/3-n) 



(1) 



where Y — (1 — f3 2 ) 1 ' 2 and h is the unit vector along the 
sight line; 



(2) 



{l-P-n) 



In Papers I & II we considered that each portion of the 
surface of the knot, denoted by the differential solid angle, 
dfl', will have an angle to the observer's line of sight that 



can range between 9 — cu/2 and 9 + cu/2. Then the values 
of 8 and v apP arising from each dQ! patch vary (unlike for 
the commonly considered case where the jet is assumed to 
be a pencil beam with 10 = 0). Here we generalize those 
results by considering the possibility that the magnitude of 
the bulk velocity, c/3, might also vary as a function of the 
angular distance (r) of that patch from the jet axis. 

The following basic relations are used for computing the 
effective values, f3 c g and 8 c g, for a knot, as would be deter- 
mined from VLBI observations which are unable to resolve 
the width of the jet (i.e., knot): 



= f s p (n')s e 

J SI 



W^Sl' = A(9)S e 



Ap{0); 



/Ipp.cff = J- / m')S p (Q')S em (Q')dn'. 

dobs J n 



(3) 



(4) 
(5) 



Here, S em and S t s are the emitted and the (Doppler 
boosted) observed flux densities, respectively, O' = /3 ■ h/ (3 
denotes the location on the knot given by the direction co- 
sine between the local velocity and the line of sight to the 
nucleus (which is no longer merely cos#), O is the entire 
solid angle subtended by the knot, A(9) is the flux boost- 
ing factor averaged over the radio knot's cross-section, and 
p is defined later in this section. Although we allow for a 
functional dependence for S em in Eqs. (3) and (5) our com- 
putations assume that S em is actually independent of Q! 
(i.e., a uniform intrinsic surface brightness across the knot). 

In principle, several peaks in a VLBI jet could belong 
to a common shock/knot yet appear separate only because 
of their different apparent velocities (owing to their different 
velocity vectors within a cone). In some cases, this approach 
has already suggested extremely large bulk Lorentz factors 
(> 100, cf., Jorstad et al. 2004); however, there is usually a 
bias against accepting and reporting such extreme speeds. 
Moreover, such associations are unlikely to be the norm, 
since there is good evidence that individual bright VLBI 
knots arise from discrete ejection events in the nucleus, since 
their trajectories can mostly be extrapolated to begin from 
the core at the times of emission of discrete flares observed 
at millimetre wavelengths (e.g., Savolainen et al. 2002). Sec- 
ondly, due to limited resolving power, the apparent motion 
of individual VLBI knots is commonly determined by mon- 
itoring the shift in the "centroid" position of a given knot 
(e.g., Cohen et al. 2007). To give expression to this practical 
situation, we have computed the flux weighted speed (Eq. 
5). Clearly, this would no longer be required when the spa- 
tial resolution and the dynamic range of the VLBI image 
undergoes a substantial improvement. 

The nominal viewing angle, 6, to the jet axis measures 
the angular offset of the circular radio disc's centre (i.e., the 
jet axis) from the direction of the AGN core (Fig. 1). The 
probability distribution of the viewing angle, V(9) (sec Eq. 
7 below), is used to compute the expectation value of 9 for 
any particular combination of T,u>,p and q, with the last 
two parameters defined below. 

The integration was carried out numerically by divid- 
ing the circular radio disc (i.e., the radio knot) of diameter 
lo and centred at 9 into small pixels of angular size equal 
to 0.01^. For each pixel (i,j) in the disc, its angular off- 
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Figure 1. Cartoon of a relativistic jet of mis-alignment angle 8 and full opening angle ui. A uniform radio knot in the jet appears as a 
disc on the plane of the sky, the centre of which marks the location of the jet's axis and therefore lies at an angular separation 8 from 
the linc-of-sight to the nucleus. 



set from the AGN core was computed and combined with 
the jet's bulk Lorentz factor at that location (see below), 
in order to find: the Doppler factor Sij, the apparent veloc- 
ity, the vector f3 a pp(i,j)\ the flux boosting factor, Aij, taken 
to be 5? .; and the product of the last two terms (which 
amounts to the apparent velocity of the pixel, weighted by 
its apparent flux). Finally, for each of these parameters, the 
average taken over all the pixels across the radio disc was 
evaluated and this was taken taken to be the effective value 
of that parameter (e.g., A e g) for the entire radio knot. Note 
that the effective value of 5 for the knot is then <5 e ff = A 1 ^ ', 
where we have usually set p = 3, widely used for compact 
discrete sources (e.g., Lind & Blandford 1985) and we have 
assumed the knot to have a flat spectral index for simplic- 
ity; taking other values for the spectral index would make 
no qualitative difference to our results. We also consider the 
case where p — 2, appropriate for smooth flows. The motion 
of each pixel is taken to be ballistic (see, e.g., Kaiser 2006). 

For the case of spine-sheath type jets, too, the (trans- 
verse) shock is assumed to cover the jet's entire cross-section 
with uniform intrinsic emissivity. In the absence of any spe- 
cific form argued for in the literature, we adopt an exponen- 
tial approximation: 

r(r) = Toe-* 9 '", (6) 

where r is the angular separation from the centre of the radio 
disc/knot (i.e., from the jet axis). We evaluated results for 
such stratified jets, taking two representative values for q 
(1 and 2). The choice of the minimum To allowed is then 
dictated by the constraint that the corresponding T(r = 
w/2) remains above the physical limit of T = 1. 

Note that the probability of finding a jet at a viewing 
angle 0, in a flux-limited sample, is given approximately by 
(Paper I; Cohen 1989): 

v(e)de qc sine aI b (6) de, (7) 

where the exponent | is the typical slope of integral source 
counts at centimetre wavelengths (e.g., Fomalont et al. 1991) 
and A e g was taken from the present computations (see 
above). Eq. (7) was numerically evaluated from — 0° to 



90° and then used to compute the expectation values (0) for 
different combinations of To, w, p and q. 

Results of the above numerical computations are dis- 
played in Figs. 2-4, for both uniform V jets (q = 0) and 
for the stratified jets (q — 1 and 2). Fig. 2 shows the ex- 
pectation values of 0, computed for a range of F, for three 
representative values of uj (1°, 5° and 10°) and for p = 2 and 
p = 3. Figs. 3 and 4 display, for p — 3 and 2, respectively, 
the dependence of /3 e g and <5 e ff on the jet opening angle uj, 
for three values of To (20, 50 and 100). Analytical expres- 
sions for the curves fitted to the computed data points are 
given inside each of the panels. Note that any given point 
along these curves refers to the corresponding expectation 
value of 9 for that point, as computed using the values of 
To, u), p and q for that point. 



3 DISCUSSION 

Fig. 2 shows a steep initial decline of {8} with increasing 
Fo, regardless of the chosen values of u or q. This decline 
continues until the regime of ultra-relativistic jets (Fo > 30) 
is reached, whereafter (9) becomes essentially independent of 
To. The near constancy of (9) for extremely relativistic jets 
is particularly striking for the conical jets having a uniform 
r and at least a moderate opening angle (a> ~ 5°-10°). 
Further, in this ultra-relativistic regime, the tendency for (9) 
to increase with uj is found to accelerate with the increase 
in the opening angle uj. Interestingly, these trends found for 
the uniform Y jets are also shared by both representative 
forms of spine-sheath jets considered here. Quantitatively, 
for uj — 5°, ultra-relativistic jets of all the three forms (i.e., 
q = 0, 1 and 2), have (8) in the narrow range from 1° to 2°. 
The corresponding range of (0) for uj = 10° jets is 2° to 4°, 
for all Ts above ~ 30. From Fig. 2, it is also evident that 
in a typical radio flux-limited sample of blazars, a larger 
opening angle of an ultra-relativistic jet would correspond 
to a considerably larger (0); hence a milder fore-shortening 
due to projection is typically expected (see Paper II for an 
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Figure 2. Expectation values of the viewing angle, {8), against the jet (central) Lorentz factor, To, for uniform (top) and transversely 
structured (middle and lower panels) jets. The left panels correspond to p = 2 and the right to p = 3. Solid, dashed and dot-dashed lines 
are for ui = 10°, 5°, and 1°, respectively. 



expanded discussion of this last point for uniform velocity 
jets). 

We now examine the influence of ui on the effective ap- 
parent speed of a radio knot, c/? e ff, and its effective Doppler 
factor, 8 e s (Eqs. 4 & 5). These are displayed in Fig. 3 (p = 3) 
and in Fig. 4 (p = 2) for the computed expectation values (8) 
(which is a function of w, To, p and S and hence varies along 
each curve, as explained above) . Good analytical fits to sim- 



ple exponential curves were found to all the computed values 
of /3 e ff and they are quoted within those figures. Exponential 
functions also gave good fits to all but two of the 18 cases 
considered for ci e ff; in those cases the displayed polynomial 
fits were found to be substantially better. 

Firstly, it is seen that the decline of /3 e s with to is 
sharper for the knots associated with jets of higher T (for 
both uniform and stratified types). On the other hand, for 
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Figure 3. The effective speeds (left) and effective Doppler factors (right) as functions of the full opening angle of the jet, with p = 3. 
The top, middle and bottom panels correspond to To = 20, 50, and 100, respectively. The X symbols, closed circles and open circles are 
the computed values for q = 0, 1,2, respectively. The solid, dashed and dot-dashed curves give the corresponding labelled exponential 
function fits to the results except that polynomial fits are given to the <5 c ff values for the To = 100 and q = 1, 2 cases. In each panel, the 
open square symbol gives the usual result for ui = and constant T. 
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Figure 4. As in Fig. 3 for p = 2, except that exponential fits are made to all results. 



well collimated jets (i.e., lo <. 0.5°), f3 c g for the uniform T 
case would typically be 1.5 to 2 times higher than for the 
q = 1 and between 2 and 4 times higher for the q — 2 (the 
two spine-sheath cases). Thus, the fastest spine component 
of the jet flow, which is near the jet axis, would be substan- 
tially concealed in the VLBI measurements. 

It is interesting to note that the much sharper fall of f3 e g 
with lo found for the uniform F case (especially in Fig. 2) 



ensures that already for modest jet opening angles, the value 
for /3 e ff of such jets would drop below the corresponding val- 
ues for both of the spine-sheath models considered here. 
For p = 3 this crossover is seen to occur near lo = 6°, 2.5° 
and 1.3°, for To = 20,50 and 100, respectively, such that 
beyond these modest jet opening angles, f3 c g rapidly ap- 
proaches mildly relativistic (or, even sub-relativistic) values 
in the case of extremely relativistic jets with uniform V. 
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This steep fall of /? e ff was first pointed out in Paper I for 
the case of uniform V jets and we find here a similar, albeit 
milder, dependence for the stratified ultra-relativistic jets as 
well (Figs. 3 and 4). It is also worth noting that while for 
well collimated ultra-relativistic jets /3 e ff is more strongly 
suppressed (relative to the uniform F case) when F is more 
peaked towards the axis (i.e., q = 2); the opposite is found 
for the jets having a significant opening angle (lo ~ 2° to 
~ 10°), with the exact cross-over value of ui depending on 
Fo and p. In summary, provided the conical jet is moder- 
ately wide (lu ~ 10°), the measured apparent speeds of the 
VLBI knots would typically remain under 10c, even if To 
were extremely large (~ 100), not only for the uniform F 
jets but even for the stratified jets. 

The right columns of Figs. 3 and 4 illustrate the in- 
dependence of the 'effective Doppler factor' (S c g) for the 
uniform F jets and for the two stratified jet forms; recall 
that <5 c ff is the cube root (for the p — 3 case), or square root 
(for p — 2) of the flux boosting factor, A c h , averaged over 
the radio knot (Eq. 4). Again, each point of these profiles 
is meant for the corresponding value of (9), as discussed 
above. It is seen that for cases of very good collimation (ui ^ 
0.5°), the uniform F jets would typically have 2 to 4 times 
larger S c s compared to the stratified jets, implying roughly 
an order-of-magnitude stronger Doppler boost. As expected 
for stratified jets, a sharper spine-sheath contrast (i.e., a 
larger q) leads to a lower 5 e g. 

Further, as for /3 c ff, a significant reduction in <5 c ff with 
u) is the typical expectation for both kinds of jet, the depen- 
dence being stronger for extremely relativistic jets, partic- 
ularly the uniform F type (Figs. 3 and 4). This leads to a 
situation where extremely relativistic jets of both uniform 
and stratified F types end up with comparable <5 c ff values, 
once to exceeds about 5°, although the S c s values for the uni- 
form jets do remain larger than those for the stratified jets 
over the wide parameter space considered here. These still 
very high values of S e g mean that rapid variability in TeV 7- 
ray emission is to be expected in any of our ultrarelativistic 
models as the variability timescale is proportional to 5~g. 
Significant changes in TeV fluxes on sub-hour timescales 
have been reported for Mrk 421 (e.g., Blazejoski et al. 2005). 

Often, VLBI observations reveal different apparent 
speeds for the radio knots in the same jet. This is usually 
interpreted by postulating that the knots reflect 'pattern 
speeds' which can be substantially different from the under- 
lying speed of the jet (e.g., Vermeulen & Cohen 1994; Cohen 
et al. 2007). In our picture, such variations can be readily un- 
derstood in terms of surface brightness distributions across 
the different knots being dissimilar. Moreover, the observed 
lack of sources with large apparent speeds but low bright- 
ness temperatures (e.g., Kovalev et al. 2005; Homan et al. 
2006) also suggests that the intrinsic pattern speed should 
be broadly correlated to the jet's bulk speed. 

As noted in Paper I, for the case of fully collimated 
jets the bulk Lorentz factor would have to quickly drop by 
1 to 2 orders of magnitude between sub-parsec (TeV) and 
parsec (radio-VLBI) scales in order to reconcile the very 
large Doppler factors inferred from the compactness argu- 
ment with the marginally superluminal (even sub-luminal) 
motion observed for the VLBI knots (Sect. 1). At present 
there is no direct observational evidence for such drastic de- 
celeration (and the concomitant dissipation) occurring on 



sub-parsec scales in blazar nuclei. Rather, for Cygnus A, 
where sub-pc radio knot motions can be measured from mil- 
limetre VLBI and the jets are close to the plane of the sky 
(so that small changes in direction would not yield signif- 
icant apparent speed changes) the evidence favors modest 
acceleration, not deceleration, on the sub-pc scale (Bach et 
al. 2005). One aim of our model is to eliminate the need for 
such a massively rapid deceleration. 

A key question is: why is the Lorentz factor dichotomy 
so striking only for TeV blazars (e.g., Piner & Edwards 
2004)? Essentially all blazars show two correlated peaks in 
their spectral energy distributions and are now usually clas- 
sified by the frequency at which the lower frequency (syn- 
chrotron) peak is strongest (e.g., Padovani & Giomni 1995; 
Sambruna, Maraschi & Urry 1996). Now, according to the 
popular scheme unifying high energy peaked (HBL) and low 
energy peaked blazars (LBL) (e.g., Sambruna et al. 1996; 
Fossati et al. 1998), TeV emission, which is an HBL char- 
acteristic, would be more common among lower luminos- 
ity blazars. This hypothesis has been supported by a recent 
study, which indicates that the synchrotron peaks for pow- 
erful blazars (the Flat-Spectrum Radio Quasars) all remain 
below the rather high energy (:> 1 keV) which is character- 
istic of HBLs and TeV blazars (Padovani 2007). 

We also recall that extragalactic jets tend to be less well 
collimated for lower luminosity sources, although this is only 
well established on the kpc scales where their opening an- 
gles can be routinely measured (e.g., Bridle 1984). Opening 
angles on sub-pc scales can be measured fairly unambigu- 
ously only for very nearby sources such as M87 (Biretta et 
al. 2002) and Centaurus A (Horiuchi et al. 2006) where the 
jets lie rather close to the plane of the sky; both of these are 
indeed relatively weak sources with wide jets. On the scale 
of nuclear jets Blandford (1993) gives a theoretical argument 
for this weak/wide correlation. Assuming it does hold, then 
the correlation between HBL properties and intrinsically 
weaker jets would mean that the jet opening angle should 
be larger for TeV emitting jets. Likely consequences of this 
are: (a) the probability of detecting TeV blazars would be 
enhanced, since, even though the TeV emission itself is prob- 
ably beamed very sharply, its effective beaming angle is in 
fact much larger, and actually more like the jet opening an- 
gle (see, Phinney 1985; Blandford 1993; Begelman, Rees & 
Sikora 1994); (b) the wider jet would mean a bigger reduc- 
tion in the apparent velocity of the radio knot (Paper I and 
Figs. 3 and 4) . We suggest that the combination of these fac- 
tors is probably responsible for the intriguing preference of 
TeV blazars to display slower motions of their VLBI knots. 
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